Vortex Structures Formed by the Interference of Sliced Condensates 
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We study the formation of vortices, vortex necklaces and vortex ring structures as a result of 
the interference of higher-dimensional Bose- Einstein condensates (BECs). This study is motivated 
by earlier theoretical results pertaining to the formation of dark solitons by interfering quasi one- 
dimensional BECs, as well as recent experiments demonstrating the formation of vortices by interfer- 
ing higher-dimensional BECs. Here, we demonstrate the genericity of the relevant scenario, but also 
highlight a number of additional possibilities emerging in higher-dimensional settings. A relevant 
example is, e.g., the formation of a "cage" of vortex rings surrounding the three-dimensional bulk 
of the condensed atoms. The effects of the relative phases of the different BEC fragments and the 
role of damping due to coupling with the thermal cloud are also discussed. Our predictions should 
be immediately tractable in currently existing experimental BEC setups. 



I. INTRODUCTION. 

Shortly after the realization of atomic Bose-Einstein 
condensates (BECs) a remarkable experiment was re- 
ported establishing BECs as coherent matter 
waves. This experiment demonstrated (apart from self- 
interference) the interference between two BECs confined 
in a trap, which was divided into two separate parts 
by means of a repulsive "hump" -shaped potential in- 
duced by a laser beam (usually called "light sheet" ) . The 
BECs were left to expand and overlap forming interfer- 
ence fringes, similar to the ones known in optics. Anal- 
ysis of this phenomenon [2j allowed for a quantitative 
understanding of some of its key features (such as the 
fringe separation) based on the mean- field theory [|, and 
in particular the time-dependent Gross-Pitaevskii (GP) 
equation; the latter, can be expressed in the following 
dimensionless form (see, e.g., [H), 



i— = --V 2 u + V(r)u + \ufu, 



(1) 



where u(r,t) is the condensate wavefunction (and 
\u(r, t)\ 2 is the atom density), while V(r) is the trapping 
potential. In fact, this potential is also time-dependent 
since it has initially a double- well form (and the conden- 
sate is allowed to relax to its ground state) and subse- 
quently, at t = 0, the "hump" separating the two wells is 
lifted; this way, the two fractions of the BEC are allowed 
to interfere and produce the beautiful experimental pic- 
tures observed (see, e.g.. the observed pattern in Fig. 2 of 
Ref. [H). Such experiments and relevant theoretical stud- 
ies are of particular value in this setting as they allow the 
study of quantum phenomena at the mesoscopic scale, 
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but with an interesting additional twist: while the under- 
lying quantum processes are purely linear, in the mean- 
field picture, inter-atomic interactions are accounted for 
through an effective nonlinearity [see the last term in 
Eq. |T])] that significantly enriches the linear behavior. 

An important modification of the linear behavior in- 
troduced by the above mentioned nonlinearity is that the 
underlying nonlinear GP model supports "fundamental" 
nonlinear structures in the form of matter-wave solitons 
and vortices. Importantly, such structures .and particu- 
larly dark bright ((J and gap solitons p}, as well as 
vortices § and vortex lattices have been observed in 
a series of experiments over the past decade. Moreover, 
experimental observation of dynamical features of these 
structures, such as the decay of dark solitons into vortex 
structures [IfJ has also been reported. 

Here, we will focus on BECs with repulsive inter- 
atomic interactions, which support stable dark solitons 
and vortex structures, in quasi one-dimensional (ID) and 
higher-dimensional settings respectively. In that direc- 
tion, and in connection with the above setting of in- 
terfering condensates, a very interesting observation was 
originally reported in Ref. In particular, it was nu- 
merically found that the collision of two quasi-ID BEC 
fragments upon release of the light sheet may lead to 
the formation of a train of dark solitons, filling the space 
originally covered by the light sheet. The phase of the 
condensate and its jumps around the soliton locations 
offered undisputed evidence that this "nonlinear inter- 
ference" leads to the formation of the relevant localized 
nonlinear structures. Subsequent work in Ref. fl2| aimed 
to clarify the regimes where the quasi-linear interference 
of the original experiments [l[ would result, versus the 
ones where the nonlinear interference of Ref. [ll[ form- 
ing dark soliton fringes, would arise. Importantly, new 
relevant experiments were very recently reported [1 31 ] , in 
which three independent BECs interfered while trapped, 
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giving rise to the formation of vortices. The latter exper- 
iment can rather naturally be characterized as a three- 
dimensional generalization of the original proposal of 
Ref. EH]. 

In the present work, we revisit the nonlinear interfer- 
ence problem of Ref. [Ill ] in its higher-dimensional ver- 
sion, relevant to two-dimensional (2D) BECs (which are 
experimentally accessible [3]), as well to fully three- 
dimensional (3D) BECs. We use appropriately crafted 
potentials to slice the condensate in two, as well as 
in four parts in 2D settings and observe their nonlin- 
ear interference in these cases. Then, we fragment the 
3D BEC into two, four and eight pieces and examine 
the results of their merging as well. The main find- 
ing of the present work is that the nonlinear interfer- 
ence of higher-dimensional BECs typically give rise to 
the higher-dimensional analogs of the ID train of dark 
solitons. Specifically, in 2D, we find nucleation of vortex- 
antivortex pairs and vortex necklaces (which have previ- 
ously been predicted to be formed in BECs as a result of 
the snaking instability of ring dark solitons [H|]), while 
in 3D we find that the vortex patterns are replaced by 
vortex ring ones. The patterns created by the collision 
of different BEC slices become more complex for larger 
numbers of slices. In examining the robustness of the rel- 
evant mechanism, we also explore the role (in the vortex 
formation process) of the relative phases of the different 
fragments, as well as those of dissipation (emulating the 
interaction of the condensate with a non-condensed atom 
fraction) and of the time for ramping down the barrier 
separating the fragments. 

Our presentation will be structured as follows. First, 
in Sec. [TT1 we present the 2D version of the problem: the 
condensate is separated into two fragments by a light 
sheet and then allowed to collide (after light sheet re- 
moval), producing vortex structures. In Sec. [Ill wc a l so 
present results pertaining to the collision of four frag- 
ments separated by a light "cross". Next, in Sec. IIII1 we 
systematically study the 3D case and interpret the cor- 
responding findings for, respectively, two, four and eight 
condensate fragments. Finally, in Sec. IV, we summa- 
rize and present our conclusions, as well as discuss some 
interesting aspects meriting future study. 



II. TWO-DIMENSIONAL CONDENSATES 

A. Collision of two in-phase fragments 

In the 2D case, we propose two different experimentally 
feasible situations. In the more standard one, i.e., the 
direct analog of the ID case of Ref. [ll|, the potential in 
Eq. ([T]) reads 

V(x,y) = VuT(x,y) + V LS (x,y), (2) 



where, for the 2D case, 

Vwr{x,y) = ±n 2 (x 2 +y 2 ), (3) 
V LS (x,y) = V aS ech(by), (4) 

where Vwr represents the harmonic trapping potential 
(with Q, being the normalized trap strength), while Vls 
is a localized repulsive potential describing the light sheet 
(with Vq and b representing the normalized intensity and 
inverse width of the laser beam). Note that, similarly 
to the experiments of Ref. [l3j |. the presence of Vht is 
necessary to guarantee nonlinear interference of the two 
condensate fractions separated by the light sheet poten- 
tial Vls (the latter pushes the BEC atoms away from the 
vicinity of the line y = 0). 

In the case under consideration, the ground state of 
the system has the form shown in the top left panel in 
Fig. [T]a, namely two identical condensate fragments sep- 
arated by the light sheet. This ground state is obtained 
by relaxation (imaginary time integration) starting from 
a Thomas-Fermi cloud with density 

\uTF{x,y)\ 2 = max{0,^- V(x,y)}, (5) 

where fi is the chemical potential that will be set to unity 
in what follows. Our proposed experiment assumes that 
the condensate, confined in the the potential of Eq. @ 
is in the ground state at t — 0, and that both fragments 
share the same phase (see below for a description when 
fragments have deferent phases). Then, at that time, 
t = 0, we "lift" (i.e., switch-off) the light sheet by setting 
Vq — > 0, and subsequently let the system evolve according 
to the GP Eq. ([I]) . Notice that this is the direct general- 
ization in 2D of the setup proposed in Ref. [llj. During 
the evolution, the two fragments originally constituting 
the condensate expand, and eventually interfere [see later 
times in the left column of Fig. [Ha)]. In order to quan- 
tify the amount of vorticity generated by the collision of 
the different fragments we monitor the vorticity u>(r, t) [r 
corresponds to [x,y) in the 2D simulations and (x, y, z) 
in the 3D simulations] . The vorticity is calculated as the 
curl of the fluid velocity, v s , namely uj{r, t) = V x v s 
(see, e.g., [Hj]), where the fluid velocity is given by: 



The typical numerical experiment showing the collision 
between two fragments with same initial phase (Acf) = 0, 
see below) shown in Fig. [TJa) is performed for = 0.1, 
Vb = 5 and 6=1. This choice corresponds, e.g., to a 
pancake sodium BEC, containing N w 300 atoms, and 
confined in a trap with frequencies Q r = 2n x 10Hz and 
Q, z = 2tt x 100Hz. On the other hand, the light sheet 
parameters may correspond, e.g., to a blue-detuned laser 
beam providing a maximum barrier energy of ks X 24 nK. 
Note that the phenomenology that will be presented be- 
low does not change for the more realistic case of smaller 
trap strengths Q, which leads to larger numbers of atoms. 
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FIG. 1: (Color online) Collision of two condensate fragments originally separated by a quasi-one-dimensional light sheet lifted 
at t = 0. Panels (a) correspond to fragments with same initial phase (A(j> = 0) while panels (b) correspond to A(j> — it. 
For each case, the first column of panels depicts the density of the condensate \u(x, y, t)\ 2 at the times indicated (brighter 
regions indicate higher densities). The result of the release of the light sheet is the interference between the two fragments, 
see |it(£ = 8) | 2 , that eventually leads to the nucleation of vortex pairs. The second column of panels depicts the corresponding 
phases (brighter regions indicate phases close to — n and +7r while dark regions correspond to zero phase). The third column 
of panels depicts the corresponding fluid vorticity (blue/red indicates positive/negative vorticity), see text. The fourth panel 
depicts the spatio-temporal evolution of the vorticity by showing a contour slice at a tenth of the maximum vorticity. The 
evolution is responsible for the nucleation and annihilation of vortex pairs resulting in the spatio-temporal vortex filaments 
shown in the panel. In all the 2D simulations, we use, for the spatial variable, a discretization of 301 x 301 sites centered around 
(x,y) = (0,0). Also, the trap parameters for all 2D simulations are: fl = 0.1, Vo = 5, and 6=1. 



Through the interference, quasi- ID "nonlinear" fringes 
are formed, i.e., dark stripes (see atom density at t = 8 
in the left panel of Fig. [T]a) resembling dark solitons (in 
direct analogy to the ID case), but now in the 2D set- 
ting. However, it is well-known that such ID stripes are 
unstable in 2D (and 3D) towards transverse modulations 
(see, e.g., Refs. [13, EH [H, H3] and references therein). 
As a result of the ensuing snaking instability, such stripes 
break up into vortex-antivortex pairs (because the origi- 
nal zero vorticity of the solution needs to be preserved). 
This is precisely what happens in our case as well. The 
bending of the stripes [see the atom density at t = 25 in 
Fig. [TJa)] is responsible for the break up into four vortex 
pairs as is shown in vorticity plot at t = 25. As the dark 
stripes continue bending and mixing, a series of vortex 
pair nucleations and annihilations occur. 

To follow this evolution, we depict in the right panel 
of Fig. [lja) a space-time plot of the vorticity. The main 
evolution of the vorticity can be summarized as follows. 
At t ~ 15 four vortices (i.e., two vortex pairs) nucleate. 
Two more quartets of vortices are nucleated at t ps 22 
and t w 35. At t w 45 the first quartet merges and is 
annihilated. The second quartet merges at t w 58. This 
process of merging and spontaneous nucleation of vortex 



pairs and quartets remains active for long times (data 
not shown here). 



B. The role of phase in two fragments collisions 

In actual BEC experiments (cf. Ref. [l3j]) the conden- 
sate is grown inside the combined harmonic trap Vht 
and light sheet potential Vls- There are two distinct 
regimes depending on the strength of the light sheet: 
if the light sheet strength is weak, the fragments have 
enough overlap and maintain a common phase, however, 
for stronger light sheets, the fragments grow essentially 
independently carrying their own independent (random) 
phases. Another possibility in the experiment is to grow 
a condensate in the harmonic trap alone and then split 
the condensate into fragments by adiabatically ramping 
up the light sheet. Depending on the degree of adiabatic- 
ity and, more importantly, on the time the fragments are 
kept separate before releasing the light sheet, the phases 
of the different fragments will evolve independently and 
assume, effectively, random phases. Here, we explore the 
effects of such phase differences between the condensate 
fragments. We therefore consider the case where just be- 
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FIG. 2: (Color online) Effects of phenomenological damping on the collision of two condensate fragments with same initial 
phase (A(/> = 0). The panels depict the same information as in Fig. [TJa) with the addition of phenomenological damping with 
(a) 7 = 0.005 and (b) 7 = 0.1. 



fore the light sheet release the phase of the top (bottom) 
fragment is 4>i (02)- Without loss of generality we fix 
4>i = and we focus on the phase difference between the 
condensates A(j> = <f>2 — <j>i = <j>2 ■ In Fig- IHb) we depict 
the collision of two fragments as in Fig. [IJa) but with a 
different initial relative phase between the condensates of 
A<p — 7r. In the simulations, initial configurations with 
arbitrary phase differences A<j) were achieved by using 
the steady state solution found for A<p = as before, 
then applying a phase shift of A(j> to the bottom frag- 
ment (i.e., for y < 0), and running again the relaxation 
scheme (imaginary time integration) until convergence. 

The results presented in Figs. QJa) and QJb) indicate 
that the detailed evolution of the vortex formation and 
merging depends on the relative phases of the conden- 
sates. In fact, further numerical results with phase dif- 
ferences between Atf> = and A(j> — tt show similar 
vortex formation and annihilation with vortex "activity" 
decreasing (both in terms of the number of vortices pro- 
duced and the persistence time of the structures) as Acf> 
was increased from to tt (results not shown here). It 
is interesting to note that the dynamics for phase dif- 
ferences different from zero loses its four-fold symmetry 
when compared to the Ac/> = case. In all cases that we 
tried, vortex pairs are nucleated at some point in the sim- 
ulation irrespective of A(f>. However, it is evident from 
the figures that the smaller the phase difference between 
the fragments, the more vortex structures are nucleated 
and the longer they live. In fact, the extreme case of a 
A</> = 7r between fragments barely produces any vortex 
pairs [see Fig. [ljb)]. This is due to the fact that the ini- 



tial condition is close to a steady-state domain- wall (dark 
soliton stripe separating out of phase domains) and thus 
there is little extra energy for the collision of the frag- 
ments. Indeed, as it can be evidenced from Fig. [Hb) 
(where A<p = tt), the condensate always maintains a con- 
figuration similar to a domain-wall with some perturba- 
tions: the density always shows a nudge at y = and 
there is a predominant phase difference of tt between the 
upper and lower half planes. 



C. Dissipation and barrier ramping-down effects 

The GP model used above relies, by construction, on 
the mean-field description of a boson gas at extremely 
low temperatures and becomes exact at T — 0. When 
the temperature is finite, but still below the critical tem- 
perature T c for BEC formation, there exists a fraction 
of atoms that is not condensed, the so-called thermal 
cloud. This thermal cloud is in fact coupled to the con- 
densed gas and its presence produces effects that are not 
accounted for by the GP equation (cf. the insightful re- 
view in Chap. 18 of Ref. [21( and references therein). 
Phenomenologically, one of the most noticeable effects of 
the presence of a sizeable thermal cloud is the introduc- 
tion of damping (to the condensed gas). The approach 
of adding phenomenological damping to emulate ther- 
mal effects was originally proposed by Pitaevskii [22j and 
applied with a position-dependent loss rate in Ref. [23| . 
In practice, a few different implementations for the phe- 
nomenological damping are viable. Here, we follow the 



approach of Refs. [H,[13| where the GP equation is mod- 
ified by the inclusion of a damping term, thus taking the 
form 

d 1 

(1 - il)-Kj. u = - 2 + V ( r > l ) u + \ u \ 2u ~ M u > ( 7 ) 

where 7 is the damping rate (the chemical potential is 
again fixed at [i = 1). In Fig. [2] we display the collision 
of two fragments as in Fig. [lja) (i.e., same initial phase, 
Acf> = 0) for two values of damping. Figures[2ja) and[2]^b) 
correspond to the cases of weak (7 = 0.005) and strong 
(7 = 0.1) damping. For the above choices, it is relevant to 
note that in Ref. [23| it was found that the value 7 = 0.03 
corresponds to a temperature of about 0.1T C . As it can 
be observed from the figure, the case of weak damping 
[7 = 0.005, see Fig. Ufa)] behaves qualitatively the same 
as the case without damping [see Fig. [Ija)] until about 
t = 40, when a noticeable reduction in the generated vor- 
ticity sets in. This effect is even more dramatic for larger 
damping [7 = 0.1, see Fig.^b)] where it can be seen that 
after t = 65 there is a complete absence of vortices. The 
effect of damping can be easily, qualitatively, understood 
in that the added phenomenological dissipation slows the 
fragments in their collisions and, therefore, less vorticity 
is produced. 

We also studied the effect of ramp-down time of the 
laser sheet separating the different condensate fragments. 
In current experiments (see, e.g, Ref. [HI) the barrier 
created by the laser sheet can be removed using a grad- 
ual ramp-down of the laser intensity. Simulations (not 
show here) reveal that slower ramping-down times have 
similar suppressing effects on vortex formation to the 
phenomenological damping cases described above. In 
fact, this is also for similar reasons, namely that longer 
ramping-down times induce a slowing down of the frag- 
ments and hence the production of less vorticity. A 
more detailed study of the effects of ramp-down times on 
the vortex nucleation is beyond the scope of the present 
manuscript and will be reported in a separate publication 
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FIG. 3: (Color online) Similar to Fig. Ola), but for an ini- 
tial condition that contains four fragments (see \u(t = 0) | 2 ), 
being the ground state of the potential of Eq. J5). Notice 
the four-fold symmetry of the interference in this case and 
the eventual formation of a complex pattern of vortex pairs. 
Same parameters as in Fig. [T]with a = b — 1. 

space, resulting in an interference pattern with four-fold 
symmetry. The oscillation and ensuing bending of the 
resulting dark stripes results in a rich evolution of vortex 
pairs as it can be evidenced by the space-time plot of the 
vorticity in the right panel of Fig. [3] Notice that for this 
initial condition, more exotic patterns, including vortex 
necklaces (cf. also Ref. [Hj]) and structures with a higher 
number of vortices are formed. However, these patterns 
are considerably less long-lived as compared to the ones 
observed in the two in-phase fragments numerical exper- 
iment described above. A more detailed analysis of the 
role of different relative phases for the different fragments 
will be reported elsewhere [26| . 



D. Collision of four in-phase fragments 

The above cases represent 2D renderings of the setting 
as originally proposed in Ref. [TT| (which, however, pos- 
sesses a number of twists particular to 2D, as illustrated 
above). Nevertheless, a more genuinely 2D installment 
of the same "experiment" can also be envisioned and is 
proposed in Fig. [3J In this case, the light-sheet potential 
is of the form 

V LS = V (sech(aa;) + scch(6y)) , (8) 

resulting in a ground state which contains four in-phase 
condensate fragments (as is shown by the initial density 
|u(i = 0)| 2 in Fig. [3]). Subsequently, as the "light cross" is 
lifted at t — 0, the fragments attempt to fill in the empty 



III. THREE-DIMENSIONAL CONDENSATES 

We now turn to the 3D analogs of the experiments 
proposed in the previous section. We start, once again, 
by a 3D version of the ID suggestion of Ref. [ll[ , with a 
harmonic trapping potential 

V aT = ^Sl 2 (x 2 + y 2 + z% (9) 

and a light sheet potential of the form 

V L s=V sech(cz). (10) 

We use c = 1, in the typical results of Fig. [U leading 
to the fragmentation of the initial condition into two 
pieces as is shown for t = 0. Once the single light sheet 
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FIG. 4: (Color online) Collision of two 3D BEC frag- 
ments. A 3D contour plot of atomic density \u(x, y, z, t)\ 2 
at max(|u| 2 )/2.5 and the vorticity u)(x,y, z,t) at max(u;)/15 
is shown as a function of (x, y, z) at different times. Positively 
(negatively) charged vortex rings are depicted in blue (red). 
The spatial discretization is 50 x 50 x 50 sites. 




FIG. 5: (Color online) Similar to Fig. [4] but for an initial 
condensate with four fragments. In this case, both, positively 
and negatively charged vortex rings, are depicted in red. Same 
parameters as in Fig. 3] with a = c = 1. 



of the potential of Eq. (TIT)]) is lifted, interference ensues 
among the different condensate fragments (see snapshots 
at t = 60, 70, 80). As the two fragments get closer, pairs 
of vortex rings are nucleated between the fragments, ft is 
interesting to observe that these vortex rings, created in 
pairs (recall that our system conserves angular momen- 
tum) around the condensate (as smoke rings that are 
created at the edge of the flow tube), may travel inward 
and "pinch" the condensate promoting its fragmentation. 
The formation of vortex rings in our system is reminis- 
cent of the experimental realization of vortex rings in- 
duced by defects in Ref. (27[. Here, again, a complex 
process involving nucleation and annihilation of vortex 
rings persists for long times. 

We now turn to the 3D analog of our second 2D nu- 
merical experiment, using a "light cross" potential 



= (sech(cz) + sech(aa;)) . 



fill 



to split the 3D condensate into four fragments initially, 
as shown at t = in Fig. [5l These four fragments start 
interfering, producing patterns with fringes possessing 
four-fold symmetry (see snapshots at t — 18,21,25,34). 
A similar scenario, though more complex, as in the pre- 
vious experiment takes place: in particular, we observe 
the formation of vortex rings around the condensate that 
promote the "pinching" of the different interference frag- 
ments. It is very interesting to note that for longer times 
(t > 100) a recurrent evolution emerges whereby a cage of 
vortex rings surrounds the bulk of the condensate atoms 
(cf. snapshots for t — 138 and t — 208) alternated by 
turbulent-like patterns with vortices (cf. snapshots for 
t = 77 and t = 249). This behavior seems to persist for 
long times (data not shown here). 

Finally, in the same spirit as followed above, one can 
attempt to produce a genuinely 3D version of this setup, 
by splitting the condensate into 8 fragments as is shown 
in Fig. [6] This can be done through a 3D light sheet (con- 
sisting of three mutually perpendicular ID light sheets): 

Vhs = (sech(ax) 4- sech(by) + sech(az)) . (12) 

As a result of releasing the light sheets, the fragments 
interfere with an eight-fold symmetry (see snapshots at 




FIG. 6: (Color online) This last 3D case is similar to the ear- 
lier ones, but for an initial condensate with eight fragments. 
In this case also, both, positively and negatively charged vor- 
tex rings, are depicted in red. Same parameters as in Fig. [4] 
with a = b — c = 1. 



t = 36, 40 in Fig. [6]). Again we find that vortex rings are 
created in pairs around the condensate and subsequently 
migrate inward promoting the "pinching" of the inter- 
ference pattern. We also find the same recurrence phe- 
nomenon, as in the previous experiment with four frag- 
ments, whereby a cage of vortex rings surrounds the bulk 
of the condensate (cf. snapshots at t = 139, 191, 209) al- 
ternating with vortex-rich, turbulent-like states (cf. snap- 
shot at t = 247). Yet again this alternating behav- 
ior seems to persist for longer times (data not shown 
here). A general observation for the experiments with 
the two and three light sheets (in comparison with the 
more straightforward realization with one such sheet) is 
that vortex rings have a shorter lifetime and are harder 
to detect due to the complicated nature of the dynam- 
ics. This observation is, in principle, also true for the 2D 
case (as can be seen by comparing the two light sheet 
case, with that of a single light sheet). Hence, perhaps 
the most robust and straightforward configuration for ob- 
serving the nonlinear interference dynamics and the for- 
mation of vortices and vortex rings respectively in 2D 
and 3D consists of the quasi- ID configuration with two 
fragments (and a single light sheet) in each case. 
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IV. CONCLUSIONS AND CHALLENGES 

In the present work, we have discussed a series of nu- 
merical experiments, representing generalizations of the 
original suggestion of Ref. pH for a "nonlinear interfer- 
ence" , leading to the formation of fundamental nonlin- 
ear structures. Importantly, relevant experiments with 
higher-dimensional BECs (in a slightly different config- 
uration than the one proposed here) demonstrating this 
concept have been reported very recently [l3j . 

While in the ID example of Ref. [HI], the ensuing non- 
linear waves were robust dark solitons, in 2D and 3D 
settings the situation is considerably more complex, even 
though fundamental nonlinear structures still emerge. In 
the 2D setting, corresponding to a pancake condensate, 
we find that slicing the condensate with a light sheet or a 
light cross gives rise to interference patterns and the nu- 
cleation of vortex pairs. This is the situation most rem- 
iniscent of the experiment of Ref. [l3j (although the lat- 
ter involved the merger of three fragments), where vortex 
creation was observed as well. In our case, subsequently, 
a complex cascade of vortex pair annihilations and nucle- 
ations was found to persist for long times. Some vortex 
"necklaces" were also seen to be stable for long times. 
This is particularly the case for a BEC fragmented in two 
pieces, whereas the vortex patterns in the four fragments 
case have shorter lifetimes, and thus it is expected that it 
would be more difficult to observe them experimentally. 

In the simpler case of two fragments in 2D, we stud- 
ied the effects of different phases between the condensate 
fragments. An interesting extension of the results pre- 
sented herein is the study of the role of the phase differ- 



ence between fragments in a condensate with more than 
two fragments [26[ . We also analyzed the suppressing ef- 
fects of phenomenological damping (induced by coupling 
with the thermal cloud) and briefly discussed the sim- 
ilar role of the ramping-down of the light sheet on the 
formation of vorticity. 

The 3D setting is particularly interesting, since the 2D 
vortex patterns are replaced by vortex ring ones. We find 
that these vortex rings, upon nucleating in pairs around 
the condensate, migrate inward promoting the "pinch- 
ing" of the different interference patterns. Subsequently 
the condensate evolves under a complex pattern of nuclc- 
ation and annihilation of vortex rings. More interesting, 
however, is the observation that for longer times the con- 
densate seems to alternate between a pattern consisting 
of a vortex ring cage around the bulk of the condensed 
atoms and a pattern of turbulent-like vorticity. 

An immediately interesting ramification of our study 
would be to compare/contrast the situation where the 
condensate fragments originate from one component to 
that where they may originate from different components 
[28| . Such interference experiments, in part, already ex- 
ist and have interesting consequences regarding the 
character of the interaction of the two-components. It 
would be interesting to examine whether such interac- 
tions can lead to (possibly anti-correlated) vortices and 
vortex rings in the two-component condensate. Such 
studies are currently in progress and will be reported 
elsewhere. 
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